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Null test of the cosmic curvature using H(z) and supernovae data 
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We introduce a model-independent approach to the null test of the cosmic curvature which is 
geometrically related to the Hubble parameter H(z) and luminosity distance dz(z). Combining the 
independent observations of H(z) and dr(z), we use the model-independent smoothing technique, 


| : i 0 
Gaussian processes, to reconstruct them and determine the cosmic curvature ag 


in the null test 


relation. The null test is totally geometrical and does not assume any cosmological model. We show 
that the cosmic curvature og = 0 is consistent with current observational data sets, falling within 
the lo limit. To demonstrate the effect on the precision of the null test, we produce a series of 


simulated data of the models with different Qu, 


Future observations in better quality can provide 


a greater improvement to constrain or refute the flat universe with og =Q: 


PACS numbers: 98.80.-k, 98.80.Es, 98.80Jk 


I. INTRODUCTION 


Whether the space of our Universe is open, flat, or 
closed is one of the most fundamental problems in mod- 
ern cosmology. The spatial topology of the Universe is 
also closely related to other important problems such as 
the evolution of the Universe, the property of the dark 
energy, etc. The effect of allowing nonzero curvature on 
dark energy models has been studied in many papers, see, 
for example, Refs. HHA. And the detection of a signifi- 


cant deviation from Q9 — 0 would have profound con- 
sequences for inflation models and fundamental physics. 
A lot of attention has been attracted to this issue 
The ACDM model is consistent with all of the data and 
a flat universe is preferred even in recent Planck 2015 
results 


However, almost all of the papers studying and con- 
straining the cosmic curvature assume some specific mod- 
els for dark energy such as the equation of state w(z). 
These are all model-dependent and indirect methods. 
Note the fact that there is some degeneracy between 
the spatial curvature and the equation of state of dark 
energy in these studies. It would be better to detect 
the space curvature of the Universe by geometrical and 
model-independent methods. In this paper, we use a 
geometrical relation among the cosmic curvature og ), 
Hubble parameter H(z), and luminosity distance d NE 
by combining the Hubble rate and luminosity distance, 
we are able to directly determine and test whether the 
cosmic curvature deviates from 0 Lo man For this 
goal, we should focus on two independent Aum 
that directly give H(z) and dr (2), respectively. For H(z) 
data, it can be derived from differential ages of galaxies 
[^cosmic chronometer (CC) "] and from the radial baryon 
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acoustic oscillation (BAO) scale in the galaxy distribu- 
tion. As for dr (z), we use the SNIa Union 2.1 data sets. 
We use the model-independent method Gaussian pro- 
cesses (GP) for smoothing the observational data. The 
advantage of these methods is that they are all model- 
independent; hence, we need not assume any models in- 
volving dark energy and gravity theory. They are purely 
geometrical and are constrained directly by observational 
data. However, the precision of the null test is also lim- 
ited by the quality of observational data. If we want to 
detect a tiny cosmic curvature more precisely, a better 
quality of the observational data sets are also required, 
which will be also discussed in this paper. 

This paper is organized as follows. In Sec. [LI] we intro- 
duce the theoretical method for the null test of the cosmic 
curvature. In Sec. we first give a brief introduction 
of Gaussian processes, and then apply GP method to the 
null test of the cosmic curvature using two independent 
data sets: CC+BAO, Union 2.1. Further it is followed 
by a series of simulated data tests. We give discussions 
and conclusions in Sec. 


II. THEORETICAL METHOD 


In a FLRW universe, the luminosity distance dy can 
be expressed as 


er 


where E(z) = H(z)/Ho, au = —Kc?/(aoHo)?, and 
K = +1,—1,0 corresponds to a closed, open, and flat 
universe, respectively. 

Differentiating Eq. (1) and writing D(z) = (Ho/c)(1 + 
2) "dr (2) as the normalized comoving distance, we can 
obtain 


o9 


(1) 


E?(z)D?(z) - 1 


AG = DX) . (2) 
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We can see that the cosmic curvature ag can be directly 
determined by using the Hubble parameter and luminos- 
ity distance from Eq. (2). Thus, the null test of Q9 is 
straightforward. Note the fact that D(0) — 0 will bring 
about a singularity at z — 0. Therefore for more succinct, 
we transform Eq. (2) to 


De) 


EG)D'G)41 = E(z)D'(z) — 1. (3) 


Since the left-hand side (Ihs) of Eq. (3) is a nonzero when 


z z 0 if Q9 is nonvanishing, the null test of cosmic 


curvature ag is equivalent to the null test of the whole 


Ihs of Eq. (3). If we define 


e DAE) 
Ox(z) = BED) +1" (4) 
then a flat universe implies 
Ox(z) = E(z)D'(z) -120 (5) 


is always true at any redshift. A deviation from it will 
indicate a signal of nonvanishing cosmic curvature. Note 
that the theoretical value of Ox(z) is always zero at z = 
0. So the null test is just to check whether there exists 
signal of nonvanishing Ox (z) at nonzero redshifts. 

Thus we should use current observational data sets to 
reconstruct E(z) and D'(z), independently, and combine 
these two reconstructions to test whether the relation 
E(z)D'(z) — 1 = 0 holds at arbitrary redshifts. We want 
to stress here that the null test is totally geometrical and 
cosmological model independent, the reconstructions of 
E(z) and D'(z) are directly derived from observational 
data sets. 


II. NULL TEST USING H(z) AND 
SUPERNOVAE DATA 


Given some observational data sets, it is crucial to 
use a model-independent method to reconstruct E(z), 
D(z), and its derivative D'(z). There are many dif- 
ferent methodologies to reconstruct functions from the 
data. For a brief analysis, see [13]. Since we should 
use a nonparametric approach to smooth the data and 
to reconstruct the derivative, the so-called Gaussian pro- 
cesses are very suitable for our purpose. 


A. Gaussian processes 


The Gaussian processes allow one to reconstruct a 
function from data without assuming a parametrization 
for it. Here, we use the Gaussian processes in Python 
(GaPP) [16]. This GP code has been used in various 
papers for different studies The distribution 
over functions provided by GP is suitable to describe the 


observational data. At each point z, the reconstructed 
function f(z) is also a Gaussian distribution with a mean 
value and Gaussian error. The functions at different 
points z and are related by a covariance function k(z, 2), 
which only depends on a set of hyperparameters / and 
of. Here, l gives a measure of the coherence length of 
the correlation in the x direction and of denotes the over- 
all amplitude of the correlation in the y direction. Both 
of them will be optimized by GP with the observational 
data sets. In contrast to actual parameters, GP does not 
specify the form of the reconstructed function. Instead, 
it characterizes the typical changes of the function. The 
detailed analysis and description of the GP method can 
be found in fig [19]. 


z 


B. Hubble rate data and Union 2.1 


Following [LA E we proceed to an analysis based 
on observational Hubble data compiled from several 
sources, independent of SNela. We combine measure- 
ments of H(z) obtained with two methods. One is cos- 
mic chronometers, which are mainly passively evolving 
galaxies. There are 21 data points compiled by Moresco 
et al. BEI Bel. The other is radial baryon acoustic oscil- 
lations from galaxy clustering in redshift surveys, which 
gives seven data points of Hubble parameters from dif- 
ferent experiments [27 39. We summarize the total 28 
data points in Table I] 

We normalize H(z) using Hy = 70km/(s Mpc); thus, 
we get the observational data points of E(z). Then we 
can use GP method to reconstruct E(z). Note that Hp is 
just a normalization factor, whose value will not influence 
our null test Eq. (5). 

To reconstruct D(z), we use SNela Union 2.1 data 
sets [31], which contain 580 SNela data. We transform 
the distance modulus m — M given in the data set to D 
using 


m — M + 5log EJ —25—5log[((1--z)D]. (6) 
For consistency, here we also use Hy = 70km/(s Mpc) to 
normalize dr (z). Obtaining these 580 observational data 
points of D(z), we can also use GP method to reconstruct 
D(z) and its derivative D'(z). Finally, we combine the 
reconstructions of E(z) and D'(z) and apply them to the 
null test of Ox(z) in Eq. (5). We stress again that the 
null test is model independent, so we need not assume 
any cosmological model, and the two data sets of H(z) 
and supernovae are also independent of each other. 


C. Null test 


Having obtained the total 28 data points of E(z) and 
580 points of D(z), we now use the GP method to recon- 
struct them, respectively. 
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FIG. 1: Gaussian precess reconstruction of E(z) (left) from CC+BAO, D(z) (middle), and D'(z) (right) from Union 2.1. The 
shaded blue regions are the 68% and 95% C.L. of the reconstruction. The flat ACDM model (red line) with Qmo = 0.27 is also 


shown. 


Index 2 H(z) Refs. 
1 0.090 69+12 [25] 
2 0170 83+8 [25] 
3 0.179 75+ [25] 
4 0.199 75 + 25] 
5 0.240 79.69 + 2.32 [27] 
6 0.270 77+14 [25] 
7 0.350 821+49 [28] 
8 0.352 83+14 [25] 
9 0400 95+17 [25] 
10 0.430 86.45 +3.27 [27] 
11 0.440 82.6+7.8 [29] 
12 0.480 97462 [25] 
13 0.570 92.4+4.5 [30] 
14 0.593 104+13 [25] 
15 0.600 87.9+6.1 [29] 
16 0.680 92+8 [25] 
17 0.730 97.347 [29] 
18 0.781 105+12 [25] 
19 0.875 125+17 [25] 
20 0.880 90+40 [25] 
21 0.900 117+23 [25] 

22 1.037 154420 [B5] 
23 1.300 168+17 [25] 
24 1.363 160+33.6 [26] 
25 1.430 177418 [25] 
26 1.530 140+14 [5] 
27 1.750 202+40 [5] 
28 1.965 186.5 +50.4 [26] 


TABLE I: H(z) measurements from different surveys using 
passively evolving galaxies and radial BAO. 


We can see from Fig. [I] that all of the reconstructions 
of E(z), D(z), and D'(z) are consistent very well with 
the flat ACDM model, which we assume for comparison. 
The dashed blue line is the mean of the reconstruction 
and the shaded blue regions are the 68% and 95% confi- 
dence level (C.L.) of the reconstruction. As expected, at 
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FIG. 2: Reconstruction of Ox(z) from CC+BAO and Union 
2.1. The shaded blue regions are the 68% and 95% C.L. of the 
reconstruction. The red line corresponds to the flat universe 
oy eg 


higher redshifts, the errors become large due to the poor 
quality data in that region. Using the reconstructions of 
E(z) and D'(z), we apply Monte Carlo sampling to de- 
termine the Ox (2) in Eq. (5) at each point z, which we 
want to reconstruct. In Fig.[2] it is shown that the recon- 
structed Ox (z) is consistent with the vanishing cosmic 
curvature, falling in the 10 limit. It tells us that there is 
no significant signal to indicate the deviation of the cos- 
mic curvature ag from 0 at the current observational 
data [H(z) and supernovae] level. In addition, let us 
mention that the mean value of the cosmic curvature is 
negative in the high redshift region, which is also consis- 
tent with the results from model-dependent constraints 
in the literature. 


D. Mock data 


To demonstrate how a large number of data with what 
accuracy of the error will affect our null test when re- 
constructing E(z) and D(z), we firstly simulate a data 
set of 128 points for E(z). Adopting the methodol- 
ogy in [23], we draw the error from a Gaussian dis- 
tribution: og ~ N(G,e) with å = (04 + o_)/2 and 
€ = (o. —0-)/4, where 04 and o_ are the two straight 
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FIG. 3: Gaussian process reconstruction of E(z), D(z) (top), D'(z), and reconstruction of Ox(z)(bottom) obtained from a 


mock data set of E(z) and future DES, assuming the concordance model with ng = 0 (red line). The dashed blue line is the 
mean of the reconstruction, and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, respectively. 


lines that bound the uncertainties o(z) of the observa- 
tional E(z) data from above and below, respectively. 
Then E(2)sim is sampled from the Gaussian distribution 
E(z) sim ~ N(E(z) fia, o E), where E(z) pia is the theoret- 
ical value from the fiducial model. 

As for simulated D(z) data, we create mock data sets 
of future SNela according to the Dark Energy Survey 
(DES) [33]. The DES is expected to obtain high qual- 
ity light curves for about 4000 SNelas from z — 0.05 
to z — 1.2. From Table 14 in [33], we can calcu- 
late the errors of D, op, and the corresponding num- 
bers of SNela for each redshift bin. At every redshift 
point z, D(z)sim is sampled from the normal distribu- 
tion D(z)sim ~ N(D(z) fia, 0 p). 

Having obtained the simulated data sets of E(z) and 
D(z), we use GP method to reconstruct E(z), D(z), and 
D'(z). Then we combine the reconstructions of E(z) 
and D'(z), using Monte Carlo sampling to determine the 
Oy (z) and check whether GP can recover its theoretical 
value and distinguish it from ag =0. 

We now simulate the data points for three different 
fiducial models: concordance model, namely, ACDM 
model with ag = 0 and Nm = 0.3; two models with non- 


vanishing cosmic curvature, ag = +0.16 and Q,, = 0.3. 
We want to check whether the GP method can detect or 
recover all of them and distinguish them from each other. 
The results are shown in Figs. BIB] respectively. 

We can see from Fig. B]that E(z), D(z), and D'(z) are 
all reconstructed very well from the mock data sets as- 


suming the concordance model ag =0. And the recon- 


structed Ox (z) is also consistent with the concordance 
model nicely. Ås expected, at higher redshifts, the errors 
become large due to the poor quality data in that region. 
Furthermore, we see from Figs. [4] and [5] that the recon- 
structions of model og = +0.16 also recover the fiducial 
model very well, falling in the 10 limit and obviously de- 
viating from the concordance model av = 0 at 95% 
C.L. The large errors and not so good reconstructions at 
high redshifts (z > 1.0) are due to the poor quality data 
in that region for there are no simulated data of D(z) at 
z > 1.2. Anyway, we can claim that with these quantities 
and errors of the observational data for E(z) and D(z), 
the GP method at least has the ability to detect the cos- 


mic curvature og > 0.16 and to rule out og = 0 at 20 
C.L. 


However, as we have pointed out that the theoretical 


value of Ox (2) is always zero at z = 0. So the theoreti- 


cal values of Ox (2) for these models with different og 


deviates from the model with og = 0 tiny at low red- 
shifts. As z increases, the difference becomes large. As a 
result, we can see from Figs. [4Jand[5] that it is very hard 
to rule out og = 0 at low redshifts. Maybe it can work 
out at a middle redshift (0.6 < z < 1.0) as Figs. [4] and [5] 


show, but it is not always helpful when ao is smaller. 


So, if we want to detect a model with og < 0.16 or even 
smaller and can rule out ag = 0 at 95% C.L., a larger 
number and higher quality data sets are required. 


Figure|6]shows the reconstructed Ox (z) also using the 
mock data sets of E(z) and D(z) but with a quarter of 
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FIG. 4: Gaussian process reconstruction of E(z), D(z) (top), D'(z), and reconstruction of Ox(z)(bottom) obtained from a 


mock data set of E(z) and future DES, assuming the fiducial model with og = 0.16 (green line). The dashed blue line is the 
mean of the reconstruction, and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, respectively. The 


concordance model o9 — 0 is also shown (red line). 
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FIG. 5: Gaussian process reconstruction of E(z), D(z) (top), D'(z), and reconstruction of Ox(z)(bottom) obtained from a 


mock data set of E(z) and future DES, assuming the fiducial model with ao = —0.16 (green line). The dashed blue line is 
the mean of the reconstruction, and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, respectively. The 


concordance model ao = 0 is also shown (red line). 
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FIG. 6: Reconstructions of Ox(z) for models with different ng , the errors are a quarter of the original errors of the mock 
data sets. (a) 2% = 0.05, (b) 2 = —0.05, (c) NY = 0.09, (a) RQ = —0.09, (e) Q = 0.16, (£) AY = —0.16, The shaded 


blue regions are the 68% and 95% C.L. for the reconstruction. 


the errors. The reconstructions are more nice. And the 
null test is more precise; it can detect the model with 
(| = 0.05 or even smaller, ruling out Q = 0 at 95% 
C.L. 


IV. DISCUSSIONS AND CONCLUSIONS 


In this paper, we have introduced a nonparametric ap- 
proach to making a null test of the cosmic curvature. Us- 
ing the Gaussian process method, we reconstructed the 
Hubble rate and distance-redshift relation [E(z), D(z), 


and D'(z), independently. In the reconstruction, we 
needed not to assume any cosmological model. By com- 
bining the reconstructions of E(z) and D'(z), we can 
determine Ox(z), which is related to the cosmic curva- 


ture ag. We have shown that ag = 0 is consistent 
with current data sets (CC+BAO, Union 2.1), falling 
within the lo limit, although the mean of the recon- 
structed curvature Ox is negative in the high redshift 
region, which is also consistent with the results from the 
model-dependent constraints in the literature. In addi- 
tion, note that the mean of the reconstructed curvature 
Ox in Fig. [Bl is positive in the high redshift region for 
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the fiducial flat ACDM model. In this sense, our result 
shown in FigP] indicates that there is a little possibility 
for a closed universe. 

To demonstrate how a large number of data sets with 
different accuracies and errors will affect our null test, 
we create mock data sets of E(z) and D(z) using the 
methodology proposed in Refs. and . We found 
that with current quality of simulated data, the GP 
method can recover and distinguish models with a differ- 


ent cosmic curvature with an order of 1071 from ag = (. 


However, if we want to detect even smaller ag, we 


should decrease the uncertainties of the data sets. Mak- 
ing the errors a quarter of the mock data, it can rule out 


og = 0 from the cosmic curvature with an order of 107? 
at 95% C.L. 

Based on a one-parameter extension to the six- 
parameter ACDM model, CMB data {ol [34 provide 


à strong constraint on ag. This constraint can be im- 


proved dramatically by adding BAO data that break the 
geometric degeneracy between og and Ho. However, 
the model assumption and priors on model parameters 
may bias estimates of the cosmic curvature. Here, we 
proposed a model-independent method for reconstruct- 
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